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Abstract 

A complete derivation, from first principles, of the reaction-rate formula for 

a generic reaction taking place in an out of equilibrium quantum-field system 

is given. It is shown that the formula involves no finite-volume correction. 

Each term of the reaction-rate formula represents a set of physical processes 

that contribute to the reaction under consideration. 
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I. INTRODUCTION 



Ultrarelativistic heavy-ion-collision experiments at the BNL Relativistic Heavy Ion Col- 
lider (RHIC) and at the CERN Large Hadron Collider (LHC) will soon start in anticipa- 
tion of producing quark-gluon plasma (QGP). Confirmation of the QGP formation is done 
through analyzing rates of various reactions taking place in a QGP. So far, the reaction- 
rate formula is derived for reactions taking place in the system in thermal and chemical 
equilibrium [|T]-^. The actual QGP is, however, not in equilibrium but is an expanding 
nonequilibrium system. 

In this paper, as a generalization of we present a first-principles derivation of the 

reaction-probability formula for reactions occurring in a nonequilibrium system. We find 
that the formula involves no finite-volume corrections. We also find from the procedure of 
derivation that different contributions to the reaction-probability formula have clear physical 
interpretation, which is summarized as "out-of-equilibrium cutting rules." 

In Sec. II, we derive from first principles the formula for the transition probability of a 
generic reaction taking place in a nonequilibrium system. In Sec. Ill, specializing to quasi- 
uniform systems near equilibrium or nonequilibrium quasistationary systems]^ we further 
deduce the formula, finding that the formula is written in terms of the closed-time-path for- 
malism of real-time thermal field theory In Sec. IV, we present a calculational procedure 
of a generic reaction-probability formula obtained in Sec. III. 

II. NONEQUILIBRIUM REACTION-PROBABILITY FORMULA 

A. Preliminaries 

The formalism presented in this paper can be applied to a broad class of theories including 
QCD (cf. the end of Sec. Ill), but, for simplicity of presentation, we take a system of self- 



Pramework for dealing with such systems is comprehensively discussed in |5|. 
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interacting, neutral scalars 0's with mass m and interaction. The system is inside 
a cube with volume V = L^. Employing the periodic boundary conditions, we label the 
single-particle basis by its momentum pk = 27rk/L, kj = 0, ±1, ±2, ■ ■ ■ , ±oo {j = 1, 2, 3). 
Physically interesting reactions are of the following generic type, 

{A} + nonequilibrium system — ^ {B} + anything . (2-1) 

Here {A} and {B} designate group of particles, which are different from 0. Examples are 
highly virtual particles, heavy particles, and particles interacting weakly with 0's. Gener- 
alization to more general process, where among {A} and/or {B} are </)'s, is straightforward 
(cf., [^). For definiteness, let us assume that {A} consists of / ^'s and {B} consists of /' 
$'s. Here $ is a heavy neutral scalar of mass M, so that $ is absent in the system. For 
simplicity of presentation, we assume a $-0 coupling to be of the form —g^(f)'^/n\ {n>2). 
The transition or reaction probability V of the process (|2.1|) is written as 

r = AT/V , (2.2a) 

•^^EE E E({^}; Wk}|5t|K};{i?}) 

{k} {rik} {m^} {n'-^} 

X {{B}; K}| S\{n^}; {A}) {{n^}\ p\{m^}) S , 

(2.2b) 

^ = EE E E(W|5t|K}) 

{k} {rtk} {mk} {"kl 

X(K}|5|{M)({MIPIW)S. (2.2c) 

Here S is the symmetry factor f^, p is the density matrix, and {{B}; {n[^}\ S \{ni^}', {A}) is 
a S-matrix element of the vacuum-theory process, 

{A} + {nk} - {B} + {<} , 

where {^k} denotes the group of 0's, which consists of the number nk of 0k (0 in a mode k ). 
In Eqs. (|2.2|) , I]{k} denotes summation over momentum/momenta of 0/0's in the final state 
{B}), and of 0/0's in the "two" initial states |{mk}; {^}}) and |{^k}; {A}). (Among 

3 



the final states {B})y is |0; {B}). This is also the case for "two" initial states.) Note 

that the perturbation series for T> starts from 1, 



V = ! + ■■■ . 

It is to be noted that {A} and {B} in (S), which we write {A, B}s, are not necessarily 
involved in one connected part of (S). This is also the case for {A, B}su We assume that, 
in W = {S'^){S), {A, B}s and {A, B}si are involved in one connected part Wc (g W). 
Then, W consists, in general, of Wc and other parts which are disconnected with Wc and 
include only 0's. Generalization to other cases is straightforward ||^ . It should be remarked 
on the form of p in Eqs. (|2.2|) . Let us recall the following two facts. On the one hand, the 
statistical ensemble is defined by the density matrix at the very initial time tj (~ — oo). On 
the other hand, in constructing perturbative framework, an adiabatic switching off of the 
interaction is required Then, p in Eqs. ( p.2| ) is a functional of the in-field (piniti,^.) 

that constitutes the basis of perturbation theory. 

As will be seen below, diagrammatic analysis shows that Af, Eq. ( |2.2bD , takes the form. 



^ = KonV, (2.3) 

where A/'con corresponds to a connected diagram and T> is as in Eq. (p.2c|) . Then, we have 

The S'-matrix element in vacuum theory is obtained through an application of the re- 
duction formula 



{{B};{<}\S\{n^};{A}) 

j=l m=l k 



"k 

E E '^(^k - ^k ; < - QN-^' 



2k ''k I / I' 

X n {^K,n') n (^^k,n) n n '/'n n n 

n'=l n=l n'=l n=l J j=l m=l 

where T is the time-ordering symbol and 



m 



(2.4) 



[ 






^ 1 








4! ^k! J 



1/2 



(2.5) 



In Eq. dJ) 5( 



; ■ ■ ■) denotes the Kronecker's (5-symbol and 



■•<l>(x), 
• ■ ■ ^fa;) 



(2.6) 



Here 



-iP]^-x 



V2EkF 

with = (Ek, pk) and Fj{x) [G*^ix)] the wave function of jth $ (g {A}) [mth $ (g {5})]. 
Z's in Eq. ( p. 61) are the wave-function renormahzation constants. It is to be noted that, in 
Eq. (|2.4| ), among rik {n[J of 0k's in the initial (final) state, («k) of 0k's are absorbed in 
(emitted from) the («k) vertices in S. Remaining rik — ik (= — i]J of 0k's are merely 
spectators, which reflects only on the statistical factor in jFj in Eq. (|3.19|) below. 

(S) in V in Eq. (|2.2c|) is given by a similar expression to Eq. ( |2.4| ), where factors related 
to the $ fields are deleted. 

/^From the form for {S), Eq. ( p.4|) , we see that the permutation of 0„ (n = 1, ■ ■ ■ , Zk) and 
the permutation of 0„' (n' = 1, ■ ■ ■ , ^k) give the same Feynman diagram (in vacuum theory), 
and then Zk'^k' same diagrams emerge. Taking this fact into account, we may write ( |2.4| ) in 
the form, 

{{B}; {K}\ S \{n^}; {A}) 



m=l 



nc?^'^!^! n / ^^a./pja,) n / ^tk./^jCk, 

k \i=i / \j=l 



x^(M,{C};{0,W) 



(2.7) 



where i]^ = — nk + and A is the truncated Green function in configuration space (in 
vacuum theory), and, e.g., {y} collectively denotes yi,y2, ■ ■ ■ ,yv- 

Among the Feynman diagrams for A, are some diagrams, in which some ^'s (g {^}) [C's 
(g {C})] coincide with x's (g {x}) and/or y's (g {y}) and/or C's (g {C}) [^'s (g {E})]. In 
such cases, A is understood to include corresponding 5-functions, e.g., 5^(^ki — Xi). 

The expression for {S^), the complex conjugate of (S*), is obtained by taking the complex 
conjugate of Eq. ( p^) or Eq. ([2^71 ), where we make the substitution (cf. Eqs. ( p.2b|) and 

(PD), 



nk mk, mk (= n^) «k ^ Jk, 



4 ^ Jk- 



This applies also to the expression for (S"^) in Eq. ( p.2c|) . 

Substitution oiW = {S'^){S) into Eq. ( p.2b| ) yields, with obvious notation. 



A/'= [n/ d%d%F^{x,)F;{x'^)]j / d'ymdS'^G*^{y^)GMm) 

n {f[J ^'a,/p.(a,)) fn/^tk,/;,(Ck, 

n / d%^f;M^)\ f n / ^^c^./pJc, 



xEEEEE 

{k} {ik} {ik} 

Jk 



X 

Vi=i " / \i=i 

x5 W({x'}, {e'}; {C}, {y'} : {y}, {C}; {O, W) S . 

Here «k = ""-k " '^k + ^k, Jk = "-k ~ "^k + jk, VV = A* A, and 

'^^ E fnC?^C?^k!4!jk!jk!) (W) I P|{mk}) 

{nk} \ k / 



(2. 



(2.9) 



B. Statistical factor S 

Here, it is convenient to introduce creation and annihilation operators, a^^^ and a-p^, 
which satisfy [apk^c^p^,] = (^k.k' and [ctpk^ctpk/] = 0. A Fock space JF is constructed on |0), 
which is defined by ap^|0) = 0. For the vector | ) (g JF) that satisfies aj^^Op^l ) = n-pj\ ) 
{n-p^ = 0, 1,2, ■■■), we use the same notation as in Eq. (|2.9|) , |{nk}), since no confusion 
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arises. A key observation here is that, using the form ( p.5|) , one can easily show that S, 
Eq. (p.9|), may be represented as 



'5=E(w| in^M In 

{"k} 



1=1 



where we write 



>i=i / \i=i 



(2.10) 



{pi, ■ ■ ■ ,Pi} = { • ■ • ,Pk, ■ ■ ■ ,Pk, ■■■} 



and then i = I]k^k- Similarly, i' = I]k^k) J = Z^kA, and j' = Y^kiL- Note that ({nk}| and 
|{mk}), in between which p is sandwiched, are as in Eqs. ( |2.2| ) and (|2.9|). 

Let us write S, for short, as 5 = (6162 ■ ■ ■ b^) {N = i + j + i' + f). Let /i, ■ ■ ■ , be a 
solution in positive integers of 



Y.h=^ (1 < m < AT) 



(2.11) 



Pick out /i 6's out of 61, 62, ■ ■ ■ , and pick out ^2 out of remaining 6's, and so on, to 
make m groups. 



(2.12) 



where 1 < ii^+i < ^h+h+i < ■ ■ ■ < iN-im+i ^ ^- In Eq. ( p.l2| ), let hi and 6^ are in between 
one set of curly brackets. Then, if / < 6; is located at the left of bi' and vice versa. We 
are now in a position to write 



S{b, ■ ■ ■ 6jv) 

N 



X Sc{bi^^^^ ■ ■ ■ bi^^^^^) ■ ■ ■ Sc{bi^^^^^-^ ■ ■ ■ bi^] 



(2.13) 
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Here, the second summation X^fg runs over all solutions in integers of Eq. (|2.11|) and the third 
summation J2gv runs over all ways of making m groups as in Eq. ( p.l2| ). From Eq. ( p.l3| ), 
Sc is determined iteratively. For example, 

Scihh) = 5(6162) - 5(61)5(62) 
^0(616263) = 5(616263) - 5,(6162)5(63) - 5,(6163)5(62) 
-5(61)5,(6263) -5(61)5(62)5(63). 

Thus, we have, with obvious notation. 



S= E EE'5c(---)'5c(---)---'5,(---)- 

m=l gr 



(2.14) 



In the case of equilibrium system, all but (apCtq) and {al^ap) vanish. From the definition of 



5,, it is not difficult to show that, for > 3, 



5,(6j,6j, . . . 6jJ = 5,( : bj-^bj^ ■■■h-) 



(2.15) 



where ': • • • :' indicates to take the normal ordering with respect to the creation and anni- 
hilation operators. 



C. Reaction-probability formula 



Now, J\f in Eq. ( |2.^ ) may be written as 

^= {jlj d%d%F^ix,)F;ix'^)^ (^ilj d%d%G]{y,)G^{yr 

X E f fl / d%Y. ^ e-^^'^M f fl / d%Y.^ 



n / dx E -1^^'^'^ n / dx E 



-.e -'^f'^^ 



x5 w({x'}, {a, {('}, {y'} ■■ {y}, {Ch {O, {^}) s . 



(2.16) 



Carrying out the integration over ^'s, ('s, ^"s, ("s and the internal spacetime vertex points, 
which are included in W, we obtain, with obvious notation. 
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M= [ilj d%d%F,{x,)F;{x'^yj l^YlJ dS,d%G]{y,)G,{y'^ 

^,E'(K72^) (e^^v^J 

x5 W({x'}, {p'}; {q'}, M : {y}, {q}; {p}, {x}) S . (2.17) 
Let us Fourier transform the wave functions Fj{x), Gj{x) 



Fj(x) = /<ir,e-'«Jl'-^-lF,fe), 

G,(x) = /*,e-«H.-.G,(r,). (2,18) 



where W = {Ej, Vj) with Ej = Jrj + M^. In Eq. (plSD , of = (X^o, X^) is the space 



point, around which $'s are locahzed and Xco is the time, around which the reaction takes 
place. In general, Fj and Gj also depend on X^. 

Substituting ( p.l8| ) into Eq. ( p.l7| ) and carrying out the integration over Xj, x'j, yj and 
y'j, we obtain 

~ ' I' 



^f=[UJ dr,dr'^F,{r,)F;{r'^)\ Jl/ ds, ds'^G;{s,)G,{s'^) 



I I' i i' VI j' j 

X v^5(E - E E q - E p) ^^(E - E r;-; E p' - E q') 

x5W({r'},{p'};{q'},{s'} : {s}, {q}; {p}, {r})e^E'(^^-^^)-S'(^^-^P]-^= S . (2.19) 

Note that, when {S) (g VT) or (5*^) consists of several disconnected parts, corresponding 
(momentum-conservation) 5-function above becomes product of several 5-functions. 

The form for "D, Eq. ( p.2c| ), is given by Eq. ( p.l6| ) or Eq. ( |2.19| ), in which factors related 



to the $ fields are deleted. 

In general, M consists of several graphically disconnected parts. As assumed in Sec. IIA, 
all $'s are included in one connected parts A/'con- Other parts, which we write "D, include 



only the constituent particles 0' of the system. Then, it is obvious that TV takes the form 
M = Mcon'D (cf. Eq. ( |2.3| )). It is also obvious that P is a contribution to T) in Eqs. (^3). 
Then, such contribution does contribute to the reaction-probability P, Eq. ( p.2a{ ), as A/'con, 
which has already been dealt with in a lower-order level. Thus, computation of A/''s, which 
consist of one connected part, is sufficient. 



III. OUT-OF-EQUILIBRIUM REACTION-PROBABILITY FORMULA 

A. Preliminaries 

In this section, we restrict our concern to quasiuniform systems near equilibrium and 
nonequilibrium quasistationary systems, which we simply refer to as out-of-equilibrium sys- 
tems. Such systems are characterized by weak dependence of the reaction probabilities 
on Xc (cf. above after Eq. ( |2.18| )). More precisely, there exists a spacetime scale L^, such 
that the reaction probabilities do not appreciably depend on Xc, when Xc is in the space- 
time region \X^ — (= |AX^|) % L'^ with X^q arbitrary spacetime point. For such 
systems, the reactions are regarded as taking place in the region |X^ — X^qI ~ Going 
to the momentum space, this means that the contribution (to the reaction probability A/") 
from the state that includes "very soft" momentum |P^| ^ should be small. More 

precisely, the contribution from the summation- region in Eq. in which at least one 

momentum (out of {pj, q^, p^, q^}) is "very soft" is negligibly small.0 

^This is the case for most practical cases, which can be seen as follows. Let T be a typical 
scale(s) of the system under consideration. In the case of thermal-equilibrium system, T is the 
temperature of the system. Due to interactions, an effective mass is induced and the vacuum-theory 
mass m turns out to the effective mass Me//(Xc). In the case of m >> \fXT ^ Meff{Xc) is not 
much different from m and, for m ^ ^/XT, a tadpole diagram induces mass of 0(\/AT), so that 
Mf,ff{Xc) = 0(\/AT). \\fXT (or even AT) is the scale that characterizes reactions. We assume 
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Let us pick out (op) from S in ( |2.14|) , which appears in TV, Eq. ( p.l6|) , in the form 

1 



E 



2EpV 



(op) e 



-iP-u) 



(3.1) 



where oj stands for or The above observation shows that the quantity (|3.1|) does not 
appreciably depend on oj^, when \uj^ — X^q| < L^. This means that (op) ~ for > 1/ U 
and = Ep ^ Then, the argument at the end of the above paragraph shows that 

the contribution to A/" that include (a) can be ignored. Same reasoning shows that the 
contribution including (a^) and/or Sc{aa ■ ■ ■ a) and/or Sc{a^a^ ■ ■ ■ a^) may also be ignored. 

Recalling the identity (|2.15|) , we pick out from Eq. (|2.14|) one i5c(aj,^ ■ ■ ■ aj^^ap^.^-^ ■ ■ ■ Op^) 
(n > 3). In A/" in Eq. (|2.16|) , this factor appears in the form 

1 \ I I j 



Pi 



n «Pi' 



X exp 



(3.2) 



,1=1 i'=j+i 

where pio = Ep {I = 1, ■ ■ ■ , n) . It is not difficult to show that among the contributions to 
Af, there are contributions, whose counterparts of Eq. ( |3.2|) , together with Eq. (p.2|), can be 
united into the form 



C{{z})=t^~'S,{: (j){z,)---(f){zn) :) . 



(3.3) 



Here 



E 



2EpV 



that this scale is much larger than the "very soft" momentum scale, « \fXT (or AT).] 

Most amplitudes, when computed in perturbation theory (to be deduced below), are insensitive to 
the region \P^\ < 0(\/AT). Then, the contribution from the region \P^\ ^ XjL^ is small, since 
the phase-space volume is small. Incidentally, in the case of equilibrium thermal QED or QCD 
(m = 0), there are some quantities that diverge at infrared limits to leading order in hard-thermal- 
loop resummation scheme [§J9|. For such cases, more elaborate analysis is required. 
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where po = Ep and ': ■ • ■ :' in Eq. ( |3.3| ) indicates to take the normal ordering. As discussed 
at the beginning of this subsection, for the system under consideration, the function ( p. 21) 
does not change appreciably in the region |AZ^| ^ {Z = Y^'i^iZi/n). This leads to an 
approximate momentum conservation for the function (|3.2|): 



E^r- E Pl^\^VL^- (3.4) 
1=1 l'=j+l 



This is also the case for C{{z}) in Eq. (PT^). The conditions under which the initial correla- 
tions may be ignored are discussed in |]TU| . In the following, we ignore the initial correlations, 



inclusion of which into the formula obtained below is straightforward. 
After all this, in S in Eq. ( p. 161) , we keep only (a"'"a)'s: 



iS = y^(al, On' ) ■ ■ ■ (al/ On' ) 
^ ^\ -p,, %ii \ p., %i , / 

m,n gr j j' n+i "^+1 

^ai. a' + (o-L eta' ))■■■( ^ai a' + (o-l,, O^a' ) 



X 

^ ^Ph ) ■ ■ ■ "P'n + l ) ^^l', ^P'n ) • ■ • (^p', ^Pil ) ' (3-5) 

where j — n = j' — m and i' — m = i — n, which leads to i + j' = j + i' . 

Referring to (PTB), we use the following set-symbols throughout in the sequel: 

Vc, = V|uVf ; V| = {a;}U{y}, vf = {x'} U M , 
Ve = VfuVf ; Vf = {Ou{C}, Vf = {e'} U {C'} , 

and Vi = VfuVf with Vf [Vf^] the set of internal- vertex points in {S) [{S^)] (g W). When 
the vertex point (^j) or C,i {Q) coincides with one of the vertex points in V| (V|^), we 
include it in Vf (Vf^). At the final stage, V$ (VgU Vi) turns out to the set of external- vertex 
(internal-vertex) points of the out-of-equilibrium amplitude ( p.l9| ) representing V. 

B. Two-point function 



From Eq. ( p.l6| ) with Eq. (|3.5|), we pick out 



p,p' \/2EpVJ2Ep/V 
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where po = Ep, p'^ = Ep,, and p G {p} U {q'}, p' G {p'} U {q}, p E {^} U {('} and 
a e {^'} U {C}. Changing p and p' to 



Ph 



(p + p')/2, P =p-p' 



we get 



-iP+ip-a) 



(3.6) 



(3.7) 



xiV(X,p+) , 

iV(X;p+) = ^ e-^(^+-^-)^»e^P-^ 
p- 

^ (^p+-p-/2 «p++P-/2) , 

where X = {p + cr)/2, = -E|p+=piVx/2| and p*^ = {E^ + EJ)/2. It is worth mentioning in 
passing that one can easily derive from Eq. ( p.7|) P ■ dxN = 0. 
Now, Eq. (|3.6| ) may be written as 

^A{p,a) = ^V'Pe-^^<''-'>-^^2n9{po) 



x^bo- 



2 ,E^ + E. 



iV(X;p) 



(3.8) 



where 

As usual, we rewrite p = |p| in terms of po by using 5+(po 
we obtain 



in Eq. (|3.8|) . In doing so 



Po 



^Ej^E. 



:iV(X;p) 



1 + 



(p- 

\p\ 



-1/2 



iV(X;po,p) 



where p = p/p. Carrying out the derivative expansion (expansion with respect to 9x^) and 
keeping up to the second-order X-derivative terms, we obtain 

x27r5+(p2-m2)iV(X;po,p), (3.9) 
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where 5+(-P^ — m^) = 6{po)5{P'^ — m?) and 



l-Is|l(Vx-(vVx)^ 



4 dm? 

X, 



fv-V ^2 



+ 



iV(X;po,p). (3.10) 



Here d /dm? acts on — m^) in Eq. ( p.9| ) and v = p/po- 

In case of the system in which translation invariance holds, (aj,aq)c oc q, Eq. ( p.7| ) tells 
us that N{pq, p) = N{po, p) is the number density of a particle with momentum p. This 
allows us to interpret N{X;po, p) as the "bare" number density of a quasiparticle with p at 
the spacetime point X^^. (For more details, see [pi][.) 

C. Construction of out-of-equilibrium propagators 

So far, for simplicity of presentation, we have dealt with real-scalar-field systems. Phys- 
ical meaning of the propagators to be deduced below can be determined more transparent 
manner by employing a complex-scalar-field systems, which we deal with in the sequel of 
this section. Let ap (aj,) be an annihilation [a creation] operator for a particle of momen- 
tum p. The antiparticle counterpart of ap (aj^) is bp (fej,). For simplicity, we assume that 
the density-matrix operator p commutes with charge operator Q, [p,Q] =0. Then, all but 
{a^pttq), (&p&q), (ctp&q), (op&q) vauish. Same reasoning as at the beginning of this section 
shows that (a^6]^) and (ctp&q) are negligibly small. Thus, we are left with (aj,aq)'s and 

A) Let us take a Feynman diagram JF for A/" (cf. Eq. ( |2.16| )), and pick out from JF 
a vacuum-theory propagator iA^'^\zi — Zq) = {0\T(l){zi)(f)'^{z2) |0) G {S){e M). Then, we 
pick up the following two diagrams for A/". The first one is the same as JF, except that 
iA(°)(zi — Z2) is replaced by 

p j2EpV q J2EqV 
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which is involved in Eq. (|2.16|) . The second one is the same as T ^ except that il\^^^{^z\ — Z2) 
is replaced by 

q J2EgV p J2EpV 



with P = {Ep, — p), etc. Adding the above two contributions to the original contribution, 
and Fourier transforming on Zi — Z2, we obtain for the relevant part, 

i 



+ 2vr5(P^ -m^)ivf ^^;po,p) • (3.11) 



Here, N with po > is as in Eq. (|3.1CI|) with (|3.71 ), while, for po < 0, iV takes the same form 



( p.lOj ) where is defined, with obvious notation, as 



E 

p- 



i{E+-E^){zw+Z2o)/'i 



\"-p+p_/2 "-P-P-/2/ 

with, as before, E± = -Eip+ipiVx/si- 

As discussed at the end of the last subsection, N{X;pQ,-p) with po > is the "bare" 
number density of a quasiparticle with momentum p at the point X'^. Similarly, A^(X; Pq, p) 
with po < is the "bare" number density of an anti-quasiparticle with momentum — p at 

B) Starting from (5^) (g Af) that includes a vacuum-theory propagator 
and proceeding as above A), we obtain 

tA22iX;P) = [tAu{X;P)]* 
—i 



p2 _ ^2 _ ,IQ- 



+2n6{P'-m')N{^^±^;Po,p 



(3.12) 
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C) Let us take a set of Feynman diagrams T\ and Ti- T\ contains (cf. Eq. ( |3.5| )) 

1 ^ 1 



E 



X ('^q.,P. + ««P.)) (3.13) 
with zi G {C'} (5*^)) and Z2 G {C} (5)). is same as jFi except that ( |3.13|) is replaced 
by 

1 .7^ . ^ 1 



E 



E 



qfc \J'^Eq^V Pj 

with zi G {(,'} (g (5'"'')) and ^2 G {.^} (g (5)). Adding the contributions from JF^ and from 
JF2, we extract the relevant part, of which the Fourier transformation on zi — Z2 is 



iA2i{X;P) =2n5{P^-m^ 



X 



%o)+iV(^^;Po,P 



(3.14) 



D) Let us take a set of Feynman diagrams J-'[ and jFg. JF^ contains 

1 ...... ^ 1 



J2 



(3.15) 



with G {^} (g (S*)) and 2:2 £ {^'} (^ ('S'^))- -^2 is same as JF{ except that (|3.15|) is replaced 
by 

1 .7^ ^ 1 



E 



E 



iPj -zi 



X (5 

+ (^-p,^-q.)) 

with zi G {C} (e (^)) and Z2 G {('} (g (5^)). Adding the contribut ions from and from 
J^'21 W6 extract the relevant part, of which the Fourier transformation on Zi — Z2 is 

iAi2(X;P) = 2n6{P^-m^) 

9{-Po) + n(^^±^;Po,p)] . (3.16) 



Above derivation of iAij [i, j = 1,2) is self explanatory for their physical meaning or 
interpretation. The physical interpretation is summarized as generalized cutting rules, which 
is a generalization of Cutkosky's cutting rules in vacuum theory. (For more details, see @].) 
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D. Closed-time- path formalism 



iAij {i, j = 1,2) obtained above are nothing but the propagators in the closed-time- 
path (CTP) formahsm of out-of-equihbrium quantum field theory |Q. The CTP formalism 
is constructed on the directed time-path C = Ci © C2 in a complex-time plane, where 
Ci = (—00 +00) and C2 = {+00 —00). A field 0(a;o,x) with xq G Ci [xq G C2] 
is denoted by 0i(xo,x) [02(a;o,x)] and is called a type-1 [type-2] field. The interaction 
Lagrangian density is of the form, 

r - _ 

'-mt — ■'-int ''-int ) 



^ = ~{<l>u^r - A-2^^^^^^^r , = 1, 2) . 



Then, the vertex factor for the "type-1 vertex" that comes from is the same as in 
vacuum theory, while the vertex factor for the "type-2 vertex" is minus the corresponding 
"type-1 vertex factor." The CTP propagators are defined by the statistical average of the 
time-path-ordered product of fields, which are written as 

iAu{x,y) = {T(f)i{x)(f)\{y))c, 

iA22ix,y) = (T(f)2ix)(f)l{y))c= [iAn{y,x)]* , 

iAi2{x,y) = (02(z/)0i(a;))c, 

tA2i{x,y) = {Mx)4iy))c, (3.17) 

where T (T) is the time-ordering (anti-time-ordering) symbol. In computing ( p. 17] ), one 
identifies 02 with (pi. Comparing Eq. ( p.l7| ) with the above deduction of Aij = 1,2), 



Eqs. ( |3.11|) , ( |3.12|) , ( |3.14|) , and ( p.l6|) , we see that x of (f)i{x) in Eq. (|3.17|) corresponds 



to a vertex-point in (S) (g W) and x of <f)2{x) corresponds to a vertex-point in (S^). The 
vertex factors in (S) (g W) are —iX for — A(0V)^/4 interaction and —ig for —g^{(j)^ </))"' / (n!)^ 
interaction. Then, the vertex factors in (S^) (g W) are, in corresponding order to the above, 
iX and ig. This is in accord with the above-mentioned vertex factors in the CTP formalism. 
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E. Reaction-probability formula 

Observation made so far shows that in Eq. (|2.16|) with Eq. (|3.5|) corresponds to an 
amphtude in the CTP formahsm of the "process," 

E $1, + E $2, - E + E $1, . (3.18) 

j=i j=i j=i j=i 

As mentioned at the end of Sec. II, only connected J\f's contribute to the reaction-probability 
V. Thus, we finally obtain 



V= (^i[Jd%d'x'^F,{x,)F*{x'^)^ 



I' 



X 



VJ=1 



E J d^^l ■ ■ ■ ^Na {{uJk - U^fc')}) , 

(3.19) 



X 

diagrams 



where J-'i is a connected amplitude in the CTP formalism which includes all $'s. In Eq. ( p.l9| ), 
we have used {uj} for collectively denoting all the (external and internal) vertex-points and 
the summation runs over diagrams. A pair of u's, Uk and Uk', in a pair of brackets ( ■ ■ ■ ) in 
jFj denotes the vertex-points that are connected by iAki{uJk{k'),^^k'(k))- 
Here some remarks are in order. 



1) As mentioned at the beginning of section, inclusion of the initial correlations ( p.2|) or 
(Of) is straightforward. 



2) Taking the infinite-volume limit V ^ oo goes as follows: 

?^(^f/'^^' 

/(27r)3 ^ ^ 
ap^W— ^a(p), etc. 

Above deduction shows that there is no finite-volume correction, in the sense that 
there do not exist extra contributions to A/", which disappear in the limit V ^ oo. It 
should be stressed that this statement holds for periodic boundary conditions. 



3) It is clear from the above deduction (cf. Subsecs. B and C) that the CTP formahsm 
here is formulated in terms of the "bare" number density of quasiparticles. A canonical 
CTP formalism is formulated in terms of the physical or observed number density of 
quasiparticles. How to translate the former into the latter is discussed in ||10||. 



Finally, we make a comment on gauge theories. If we choose a physical gauge like the 
Coulomb gauge or the Landshoff-Rebhan variant |]ll| of a covariant gauge, the gauge boson 
may be dealt with in a similar manner to the above scalar-field case. If we adopt a traditional 
covariant gauge, a straightforward modification is necessary. 

IV. COMPUTATIONAL PROCEDURE 

In this section, we present a concrete procedure of computing the reaction probability V 
up to rath-order terms with respect to the derivatives. 
1) From jFj in Eq. ( |3.19| ), we pick out iAij{p,a), 

Since jFj includes <l>'s, the vertex-point p [a] is connected^ with a vertex-point v [v'] G V$ 
(cf. Fig. 1): 

p + a 1 



2 



i=o 
k' 

i=o 



(4.2) 



where Uq = p, = a, ujk+i = uj'i^i^i = v', with v, v' G V$. In Eq. (^I^), each pair of 
spacetime points in a pair of brackets, tUj+i and Uj [u'^ and ^+i], is connected by one or 



^Note that, in general, the vertex-points v and v' are not uniquely singled out. {v can coincides 
with v' .) However, different choices of v and v' leads to the same reaction probability V within the 
accuracy under consideration. 
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several iAki{Ljj+i,ujj) [iAk'i'{ujj,Ljj_^_{)] in jFj (cf. Fig. 1). Here, we note that v and v' may 
be written as 



V = X + v, v' = X + v' 



(4.3) 



where X is the mid-point of the external- vertex points, around which the reaction is taking 
place: 

1 



X 



2{l + l') 



I I' 



2) Using Eqs. §3) and we expand Aij{{p + cr)/2; P) in Eq. as 



1 k 

A,,(X;P) + -[-5:(^,+i-^,) 

dxA,,{X;P) + 



j=0 



k' 



j=0 



(4.4) 



where ' ■ ■ ■ ' stands for terms with higher-order derivative with respect to X. The series ( [1.4|) 
is truncated at the nth-order terms with respect to the derivatives. The approximation 
in which ' ■ ■ ■ is ignored is called the gradient approximation. 



3) Let us deal with the term with (tu^+i — Uj) in Eq. (^]5). It can easily be shown that 
{uj+i — ujj)iAki{ujj^i,ujj) becomes^ 



idP' 



kl 



UJj + UJj 



i+1 . 



P + P' 



Other terms and higher X^-derivative terms ' ■ ■ ■ ' in Eq. (|4.4| ) may be dealt with similarly. 



All other parts of JFj, Eq. (3.19), than the one (|4.1|) may be dealt with similarly. 



^As in the case of some self-energy-type subdiagram, there are several iAf^i{ujj^i,ujj)^s 
[iAf^iii{LL!j,u;j_^_-^^)^s] (cf. Fig. 1). In such a case, one chooses any one of them. 
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4) Carrying out the integrations over all vertex-points except those in V^, we have 
momentum-conservation 5-functions at each internal vertex point. 

As discussed at the beginning of Sec. Ill, the wave functions of <l>'s should be localized 
within the space region ^ U {i = 1,2,3). However, for simplicity, we assume in the sequel 
that the wave functions of $'s are of plane-wave form,^ 

G^{y) = e-'^^ V (2V^r-' + M^y^ . (4.5) 

5) We carry out the integrations over all vertex-points in V<!> to yield momentum- 
conservation 5-functions at those vertex points and we are left with integrations over the 
independent or loop momenta. Keeping the terms up to the nth-order terms with respect 
to the Xfj, derivatives, we obtain the final formula, which may be written in the form, 

V = Jd^X A{X; R[, ■ ■ ■ , R[r, ■ ■ ■ , Ri) . (4.6) 

Note that A depends weakly on X through N{X; Qk)'s. From Eq. (|4.6| ), we see that A is the 
reaction rate per unit volume. Incidentally, were it not for this X-dependence, integration 
over X in Eq. ( |4.6|) would yield VT, where V is the volume of the system and T = tf — ti 
is the time interval during which the reaction takes place. In the limit V, T — > oo, the VT 
becomes 

lim VT = (2n)^6\0). 
Example 

Here, for the purpose of illustration, we deal with the heavy-$ production process, 

^It is to be noted that, if we use the the plane- wave form ( |4.5| ) in Eq. ( |2.19| ), X-dependence 
disappears. In the procedure presented here, X-dependence of M is already (partially) taken into 
account before arriving at 4). 

21 



Out-of-equilibrium system — > $ + anything . 



(4.7) 



The system is composed of real scalar with £int = — A</)^/3!, and $ interacts with (j) 
through Ctf)<s, = —g^(f)'^/2. We analyze the contribution from Fig. 2 for V in Eqs. (^3). 
Using Eq. ( p.l6|) , we have 



M = g^X^ / rfVG'*(x') / dS'G{y') 



E 



pi y'^Ep^V 



P2 y ^-^P2 



E 



1 



V 



E 



q' \ 2EqiV 



where A*^''-' is the vacuum-theory propagator of (j) and S (cf. Eq. (2.10)) takes the form 



= ("i'l^k^^i''^ + aqaq')apiap2) • 

We compute the contributions that include only two-point functions. If necessary, the con- 
tributions including initial correlations may be written down in a straightforward manner. 
Keeping the terms that do not vanish kinematically, we have 

S = Si+S2, 

51 = (aj,,ap,)(aj,,ap2) [5q,q' + (a^aq')] , 

52 = Si 

Pl ^ P2 

We compute the contribution A/i from Si. The contribution from ^2 may be computed 
similarly. Following the procedure presented above, we obtain 

4f/ 



AAi = g^X^J d'x'G*{x') J d^y'G{y') J d^i j d^i 
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-iQ'-(y'-r). 



— z 



Here we observe that 

—2 — = -2^^^~''^^^^-y^^^^[-dP'-M') ' 

where the partial derivatives are understood to act on the "propagators" in momentum 
representation. 

Making the plane- wave approximation for G{x), 

^—iR-x 



G{x) 



v/2E$F 

we finally obtain, within the gradient approximation, 



{E^ = + M2) 



V^P, 



2TT6+{Pf -m')N{X;Pi 



X 



X 



27r6+{P^ - m^) N^X^; P2) 2tx5+{Q^ - m'){l + N{X^- Q)] 



^-^[9x2 + dxi 



1 



d 



d 



dP' dQ' 



1 



p'2 _ ^2 + ^0+ 



Q'2 -m?- zO+ 



(4.8) 



Xi = X2 = X, Q' = P' 

where X = {x' + y')/2 and P' = = Pi - P and Q = P2 + Pi - P- 

Eq. ([4.8|) corresponds to a contribution to the amplitude in the CTP formalism of the 
"process" (cf. Eq. (|3.18|) ), $2(P) ^liR), and constitutes a part of the diagram as depicted 
in Fig. 3 in the CTP formalism. As a matter of fact, Eq. ( [4.8[ ) represents Fig. 3 with 
(pio > 0,p2o > 0,qo> 0) plus Fig. 3 with (pio > 0,p2o < 0, go < 0). 



23 



REFERENCES 



[1] H. A. Weldon, Phys. Rev. D 28, 2007 (1983). 

[2] A. Niegawa, Phys. Lett. B 247, 351 (1990); N. Ashida, H. Nakkagawa, A. Niegawa, 
and H. Yokota, Ann. Phys. (N.Y.) 215, 315 (1992); 230, 161(E) (1994); Phys. Rev. D 

45, 2066 (1992); N. Asida, Int. J. Mod. Phys. A 8, 1729 (1993); tbtd., 3829(E) (1993); 
A. Niegawa and K. Takashiba, Ann. Phys. (N.Y.) 226, 293 (1993); 230, 162(E) (1994). 

[3] M. Jacob and P. V. Landshoff, Phys. Lett. B 281, 114 (1992); P. V, Landshoff, Phys. 
Lett. B 386, 291 (1996). 

[4] A. Niegawa, Phys. Rev. D 57, 1379 (1998). 

[5] K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

[6] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields (Mc Graw-Hill, New York, 
1965). 

[7] A. J. Niemi and G. W. Semenoff, Ann. Phys. (N.Y.) 152, 105 (1984). 

[8] R. D. Pisarski, Phys. Rev. Lett. 63, 1129 (1989); E. Braaten and R. D. Pisarski, NucL 
Phys. B337, 569 (1990); B339, 310 (1990); J. Prenkel and J. C. Taylor, ibid., B334, 199 
(1990). For an out-of-equihbrium generahzation, see, e.g., M. E. Carrington, H. Defu, 
and M. H. Thoma, Eur. Phys. J. C7, 347 (1999). 

[9] M. Le Bellac, Thermal Field Theory (Cambridge Univ. Press, Cambridge, 1996). 

[10] A. Niegawa, Prog. Theor. Phys. 102, 1 (1999). 

[11] P. V. Landshoff and A. Rebhan, NucL Phys. B383, 607 (1992). 



24 



FIGURES 

FIG. 1. A diagram for Ti in Eq. ( |3.19| ). i, j, k , and / are the vertex-type. Each $ is either 
type-1 or type-2. 

FIG. 2. A diagram representing AA, Eq. (|2.2b| ), for the process (4/7). The spacetime points 
^ and x' and y') are connected by a vacuum-theory propagator. The dot-dashed hne stands 
for the final-state-cut line. The group of particles on top of the figure represents the spectator 
particles. 

FIG. 3. An amplitude for the "process" ^2(-R) ~^ *I*i(^) in the OTP formalism, a part of which 
represents the contribution (^). 
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